Excitation of Spin Waves in Superconducting Ferromagnets 
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This Letter presents a theoretical analysis of propagation of spin waves in a superconducting 
ferromagnet. The surface impedance was calculated for the case when the magnetization is normal to 
the sample surface. We found the frequencies at which the impedance and the power absorption have 
singularities related to the spin wave propagation, and determined the form of these singularities. 
With a suitable choice of parameters, there is a frequency interval in which two propagating spin 
waves of the same circular polarization are generated, one of them having a negative group velocity. 
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Coexistence of superconductivity (SC) and ferromag- 
netism (FM) has been a long-standing problem since the 
beginning of the modern theory of SC J^^,^- There 
has been a renewed interest to this problem, and SC-FM 
coexistence was discovered recently in various materials 
such as ruthenocuprates 0, 0] , ZrZn2 0] , UGe2 j and 
URhGe . It was also suggested theoretically that in p- 
wave superconductors, time reversal symmetry is broken, 
and a non-zero magnetic moment arises j9j. 

The existence of the FM order in a superconductor is 
hidden by Meissner currents, which create a magnetic 
moment opposite to the spontaneous magnetic moment. 
Moreover, in the Meissner state SC eliminates another 
consequence of the FM order: the equilibrium domain 
structure However, the Meissner currents cannot 

shade such an important manifestation of the FM or- 
der as spin waves. A conventional technique to probe 
spin modes in ferromagnets, both insulators and metals, 
has been observation of Ferromagnetic Resonance (FMR) 
\n\ . There are also reports on FMR observations in ma- 
terials with SC-FM coexistence [l^ . 

In this work we theoretically investigated propagation 
of spin waves in a material, which possesses both SC and 
FM order and is irradiated by an electromagnetic (EM) 
wave. The first step in this direction was done by Ng 
and Varma , who studied the spectrum of spin waves 
interacting with vortex modes in the spontaneous vortex 
phase (mixed state in zero external magnetic field). Here 
we analyze the boundary problem for the Meissner state: 
how an incident EM wave can generate spin waves inside 
a sample. We solved the Landau-Lifshitz equation for 
the magnetization and the equations of London electro- 
dynamics, assuming that the equilibrium magnetization 
is normal to the surface and there is no external static 
magnetic field and, as a result of it, no Meissner currents 
at equilibrium. We found the spectrum of spin waves for 
two cases, depending on the stiffness of the spin system. 
If the stiffness is large enough, the spectrum looks similar 
to that of a FM insulator, with the minimum wave fre- 
quency (threshold for spin wave propagation) equal to the 
frequency of the uniform FMR, i.e. corresponding to zero 
wave vector ^| . But the opposite case of small spin stiff- 




FIG. 1: The setup geometry. 



ness is more interesting. Then propagation of spin waves 
becomes possible at frequencies lower than that of the 
uniform FMR, and waves near the threshold have finite 
wave vectors. Moreover, a monochromatic incident EM 
wave generates two propagating spin waves with slightly 
different wave vectors. One of them has a negative group 
velocity, i.e., its direction is opposite to that of the phase 
velocity. We found the surface impedance, which is in- 
timately connected with the spin wave spectrum. The 
analysis was focused on the case without explicit dissipa- 
tion terms in the dynamical equations, when all energy 
absorbed from the incident radiation is carried away by 
spin waves propagating deep into the sample bulk. But 
we discuss also the role of various dissipation sources. 
The analysis demonstrates that the microwave probing 
of spin modes, which has been so fruitful for normal fer- 
romagnets, should be also an effective technique for in- 
vestigation of unusual properties of spin modes in super- 
conducting ferromagnets (SCFM's). 

We consider a SCFM occupying the semispace z > 
as shown in Fig. ^ This is a good approximation for a 
slab so thick that one can neglect a possible reflection 
of waves from the second boundary. Uniting the micro- 
magnetism approach for ferromagnets and the London 
approximation for superconductors leads to the follow- 
ing free energy functional for our material: 

F= [ d^xl 27ra [mI + ll{diMj){diMj)] 
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where (p is the phase of the SC order parameter, $o = 
hc/2e is the magnetic flux quantum, M is the FM mag- 
netization, and M± is the magnetization component per- 
pendicular to the easy axis z, which is normal to the 
sample surface. At equilibrium the magnetization is 
Mo = Mqz {M± = 0), and the large anisotropy pa- 
rameter a > 1 (the ratio of the anisotropy energy to the 
magnetostatic energy) ensures that our configuration is 
stable against the flip of Mq into the plane parallel to the 
surface. In our geometry the static magnetic induction 
Bq vanishes both inside and outside the sample, but the 
magnetization creates a magnetic field inside the sam- 
ple: Ho — — 47rMo. In accordance with the micromag- 
netism approach, the absolute value of M is assumed to 
be constant M — AIq, and hence terms dependent on it 
were omitted. The length Id, which is of the order of 
the domain wall thickness in normal ferromagnets, char- 
acterizes the stiffness of the spin system. We consider a 
single-domain sample, since in SCFM's in the Meissner 
state there are no domains at equilibrium |lO|. 

The spin dynamics is governed by the Landau-Lifshitz 
equation 0| 
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where g is the gyromagnetic factor. Since we are con- 
cerned here with motion near the equilibrium, we can 
decompose 



M = Mo + m. 



(3) 



where m _L Mo is the dynamic part. For plane waves 



oc e 



ikz—iujt 



Eq. (0) yields: 

-5M0 X {47ra(l-f Z^fc2)m-b}, (4) 



where, since the equilibrium magnetic induction van- 
ishes, there is only a dynamic part of the induction: 
B = b. To find b, we minimize the free energy with 
respect to the vector potential A. This yields the Lon- 
don equation: 



(1 + fc^A^) b = 47rfc2A2m 



(5) 



Substituting this into Eq. I^J, we obtain the equation of 
motion for the magnetization: 



-iujm= — 47r5Mo x m 



1 + fc2A2 



with the dispersion relation 
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where ujfm = ^TragM^ is the frequency of the uniform 
FMR at zero internal magnetic field for normal FM's. In 
the limit \ ^ 00 the spin wave spectrum should trans- 
form to that of an insulator. Neglecting the spin stiffness 
[Id ^ 0) the latter can be presented as w = g{Hn + H) , 
where Ha = inMoa is the anisotropy field This 
agrees with Eq. ((TJ (at finite k) bearing in mind that 
H — —A-kMq in our geometry (see Fig. However, 
whatever large A could be, at very small fc <C 1/A the 
dynamic induction is screened out, as evident from Eq. 

so at these scales the spin- wave spectra in supercon- 
ducting and in insulating ferromagnets are different. 

The Maxwell equation dh/dt — — cV x e and the Lon- 
don equation, Eq. (jSJ, relate the EM field inside the 
sample to the magnetization: 



h = b — 47rm 



e = — b X z 

kc 



in 
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The two signs in Eq. {Tj) correspond to two senses of 
circular polarization: m^^' = m^^)(x=p«y). Only pos- 
itively polarized waves m'-"'"-', which correspond to the 
upper sign, can propagate inside the sample, and in the 
following we shall focus on this polarization. The form 
of the spectrum depends on the ratio of the domain wall 
thickness Id to the London penetration depth A, as shown 
in Fig. 121 If ^/ald > A, the minimal frequency (threshold 
for spin- wave propagation) is lo = u; fm id k — 0. On the 
other hand, if ^/ald < A, then the minimal frequency is at 
a finite wave vector, km = Xld~ 1/A2)1/2, and has 

a lower value, ujm = i^/m[l — (l/-\/a~ ^d/A)^]. This dip in 
the spin-wave spectrum was revealed by Ng and Varma 
[l3| for the spontaneous vortex phase. Waves with wave 
vectors satisfying |A:| < km have a negative group veloc- 
ity duj{k)/dk, i.e., its sign is opposite to that of the phase 
velocity uj{k)/k. Thus in this regime a positively polar- 
ized incident EM wave with the frequency between ojm 
and ujfm excites in the material two propagating waves, 
one with a positive and another with a negative group 
velocity. This unusual property is unique to SCFM's. 

Now we turn to calculation of the surface impedance, 
which specifies the boundary conditions for the EM field 
outside the sample. The dispersion relation is a quadratic 
equation for k'^ , so for a given frequency there should be 
two solutions for k'^: kf and These two solutions cor- 
respond to two modes generated by the incident EM wave 
inside the sample. Thus the situation is analogous to that 
in two mode electrodynamics, the theory suggested for 
calculation of the surface impedance of a superconductor 
in the mixed state taking into account the elastic degree 
of freedom of the vortex array 0| . In the present case 
the second mode (the second solution for k'^) originates 
from the spin degree of freedom. Moreover, in contrast 
to Ref. [iJI, where both modes inside the sample were 
evanescent (fc2 2 < 0), now dX to > ujfm one mode, and at 
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FIG. 2: Form of the spectrum in two regimes: for ^/ald/^ = 
1.03, a = 7 (solid line); and for y/ald/\ = 0.32, a = 1.01 
(dashed line). 



ujm < w < ujfm both modes are propagating. 

In order to find a proper superposition of the two 
modes inside the sample, m(^)(z) 



, we need two boundary conditions. One of 



m 

them is the usual continuity of the EM field across 
the sample boundary. Polarization of the excited spin 
wave is defined by the polarization of the incident EM 
wave outside the sample, and the surface impedance 
tensor is diagonal in the circular wave basis. In this 
basis h ~ h'^^\z) (ic =F iy) and e — e^^\z) (x ^ iy), 



where h^^^ [z 



e^^^(z)). Then the surface impedance is 



(and similarly for 
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The second boundary condition should be imposed on 
the magnetization. The simplest possible boundary con- 
dition [ll| is dm/dz = at z = 0, which means absence 
of spin currents through the sample surface. This gives 



(10) 



Together with Eq. 
the surface: 



this yields the field amplitudes at 



e(±)(0) = T« 
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1 + k^X^ k-2 
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Substituting this into Eq. we obtain after simplifi- 
cations: 



c 



(±) 



ujX'^ 



kik2{ki + ^2) 



l + {ki + ki + kik2)X^ 
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After finding the two roots kf and fc| of Eq. (|7|) we need 
to choose the signs of ki = i\/kl and k2 — ±\/Aif . If the 
wave vectors have imaginary parts, they should be posi- 
tive in order to ensure attenuation of the wave oc e*'^^~"^* 
inside the sample. But if k is purely real (propagating 



wave) , the the sign is fixed by the condition that the en- 
ergy must be carried away from the surface. One can 
check that the energy flux for each mode normal to the 
sample surface is proportional to the group velocity: 



„ . ,9 . dm c ^ 2TTauj , 
^ -A-Kalim—- + — E x H = |m 

OZ AtT l^fm 



Hence in order to carry energy away from the boundary, 
waves with a negative group velocity should have wave 
vectors directed toward the boundary. Bearing this in 
mind, we obtain from Eq. 112|) the general expression for 
the surface impedance: 



2/dA-igi/2_Q) 



IdX-^Q^/^ - 
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where Q = 1 =p Lo/ojfm, and the signs of the square roots 
are determined by the above requirements on the wave 
vectors (this is consistent with the condition ReC > 0). 

Now we shall analyze the most interesting particular 
cases, when the frequency is close to thresholds at which 
spin-wave propagation becomes possible. Two regimes 
should be considered, depending on the value of Idl X. 

1. Stiff spin system: ^/ald > X. Then the threshold 
frequency for spin- wave propagation is lu , and near the 
threshold fci ^ \k2\: 



kl 



an - A2 ' 



Since y/ald > A, the short-wavelength mode is evanescent 
(fc| < 0), while the long- wavelength mode is evanescent 
{k\ < Q) dX u) < LOfm and propagating {k\ > 0) at a; ^ 
ujfm ■ Finally the surface impedance is 



c 



(+) ^ 



kik. 



1 + fc2A2 



a{alj - A2) 

(16) 
Thus, C^+^ 

is purely imaginary below the threshold and 
purely real above it, where it grows as ^/uj — lu fm (actu- 
ally, there is also nonzero Im^^^^ above the threshold, 
but it is of higher order in small — tOfm) ■ 

2. Soft spin system: y/ald < X. In this regime the 
threshold frequency has a lower value < LOfm- Near 
the threshold \ki\ < km and |/c2| > km are close to each 
other: 



U! — UJ„ 



1 



(17) 



Also, one should take into account that the mode with 
*\ kfi^ has a negative group velocity above the thresh- 
old, so the negative sign of fci should be chosen. Then 
near the threshold is given by 



c 



(+) ^ UJmX'^ fcm(fc2 - _ 



. w„, A^ 
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Again the surface impedance is real above the threshold 
{uj > LOm) and imaginary below it (w < ujm)- 

In addition, there is still a singularity at uj — ujfm- We 
can use Eqs. l(T31) and but now, since \/ald < A, the 
short- wavelength mode became propagating (/c| > 0). 
The long-wavelength mode is now propagating {kf > 0) 
at w < LOfm and evanescent (fcj < 0) at > w/m- Thus 
at w — LOfm there is a transition from two propagating 
modes to one. On the other hand, at > w/™ Eq. flfil 
yields a purely imaginary impedance despite the presence 
of a propagating mode with real ^2. In order to obtain 
the real part of the surface impedance we need to go to 
next order of the expansion in the small parameter fci /fc2 : 



ReC^+Hc^>^/™)^- 



c(l + fc2A2)2 



Xcy/X^-al'j 



Hence at a; — uifm, 
~ y'tofm — UJ below ujf' 



to 



(19) 

changes its behavior from 
^ {uj - ujfm) above ujfm- 
Our analysis can be easily generalized to take into ac- 
count dissipation due to normal currents. For this, a 
normal current should be added to the SC current in 
the Maxwell equations. As a result, a renormalized com- 
plex A should be used instead of the London penetration 
depth: A^^ = ^-2 _ 2iS~^, where S — cj^f2i^a^ is 
the skin depth, and (T„ is the normal conductivity. Ex- 
panding with respect to the small ratio A^/J^, we obtain 
instead of Eq. lfTB|l : 



c 



(+) 



n 



1 



(20) 

In the opposite limit A — > oo, A^ — and we ob- 

tain the surface impedance for normal FM metals |l5l |. 
However, unlike in SCFM's, in FM metals the dissipation 
is an essential part of spin dynamics, which can not be 
assumed small. 

Note that at uj — u)f„i [see Eq. (|16ll ]. vanishes even 
after renormalization of the penetration depth, since at 
uniform FMR no currents are generated, so the dissipa- 
tion due to normal currents is ineffective. Then, other 
sources of dissipation become important, in particular, 
the transverse spin relaxation, which is well known from 
studies of normal ferromagnets llj. 

Our phenomenological approach is valid for any mate- 
rial with coexisting SC and FM order parameters, in- 
dependent of their microscopic origins. Moreover, we 
even believe that it is not essential whether FM originates 
from electron spins or from the orbital moment of Cooper 
pairs, which characterizes the p-wave pairing in some ma- 
terials 01 . In the latter case one cannot call the modes 
which were analyzed here, "spin waves" : they should be 
called "orbital waves" as in ^He physics. It is known from 
•^He studies that the dynamics of the orbital moment is 
similar to spin dynamics [l^ . In the case of electron pair- 



ing, when both the spin and orbital motion contribute to 
the magnetic moment, one can consider the dynamics 
of the total magnetic moment. Consequently, the cor- 
responding excitations will be "magnetization waves". 
In the case of coexistence of SC with orbital FM the 
term "ferromagnetism" is sometimes avoided, while the 
corresponding category of materials is commonly called 
"superconductors with broken time-reversal symmetry" . 
However, independently of this semantic variance, our 
analysis is relevant for these materials. 

In conclusion, we have demonstrated that propagating 
spin waves are possible in superconducting ferromagnets 
in the Meissner state, similar to those in insulating ferro- 
magnets. In contrast to the latter, in SCFM's with low 
spin stiffness there is a frequency interval, in which there 
are two propagating spin modes of the same circular po- 
larization, one of them having a negative group velocity, 
i.e., with the direction opposite to that of the phase ve- 
locity. We solved the boundary problem and calculated 
the surface impedance for the case of the easy magnetiza- 
tion axis normal to the surface of the sample. Revealed 
singularities of the dependence of the surface impedance 
on the frequency provide experimentalists with a tool for 
probing spin- wave dynamics by microwave measurements 
of materials with coexisting FM and SC. 
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